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Abstract 

We numerically analyse the evolution of the flavor non-singlet gi struc¬ 
ture function taking into account the all-order resummation of a* In^ x terms 
which is expected to have much stronger effects than the DGLAP evolution 
in the small x region. We include a part of the next-to-leading logarithmic 
corrections coming from the resummed “coefficient function” which are not 
considered in the calculation of Bliimlein and Vogt to respect the factoriza¬ 
tion scheme independence. It is pointed out that the resummed coefficient 
function gives unexpectedly large suppression factor over the experimen¬ 
tally accessible range of x and Q^. This fact implies that the next-to-leading 
logarithmic contributions are very important for the gi structure function. 
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1 Introduction 


Recent new data for the unpolarized Deep Inelastic Scattering from HERA provide us 
with much information on the quark-gluon structure of Nucleon. The HERA experi¬ 
ments cover much broader kinematical regions than before. Especially, the behavior of 
the structure function at small values of the Bjorken variable x receives much attention 
of the physicists [|I|. The small x region corresponds to the Regge limit. So we naively 
expect that the soft physics (Regge theory) may explain the small x behavior of the 
structure function. However the steep rise of the structure function in this region ob¬ 
served by the HERA experiments contradicts with this naive expectation. The physics 
at small x is now one of the most interesting subjects and many people believe that 
this problem could be handled in the context of the QCD perturbation theory [|^ . The 
approach based on the Balitskii-Fadin-Kuraev-Lipatov (BFKL) [@] equation or on the 
high-energy factorization P] seems to be very promising and it is important to reveal 
the role of the so-called BFKL Pomeron. 

In the case of the polarized structure function gi^ we have not yet had data at very 
small X. However the recent data show some rise of gi in the small x region Q]. This 
behavior again seems to contradict with naive Regge prediction (^fi ~ 0 < a < 

0.5) 0. In fact, to explain the rise of gi at small x in the framework of the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) approach, it is required to choose a steep 
function as an input parton density as far as the evolution starts at the order of 
1 GeV"^ because the evolution effect in the DGLAP equation does not produce enough 
enhancement in the small x region [^. So it is interesting to see whether the all-order 
resummation of In x terms which appear in the perturbative calculations reconcile the 
experimental behavior of gi with the naive Regge prediction. 

Some time ago, Kirschner and Lipatov ||^ considered the all order resummation of 
ttg In^ X series in the case of quark-quark forward scattering process. Recently Bartels 
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et al. have given the resummed expression for the gi structure function by using 
the Infra-Red Evolution Equation. They claim that the resummation effects may lead 
to 10 times larger results than the DGLAP ones. This analysis suggests the above 
possibility that the small x behavior of gi is explained naturally by combining a flat 
input (non-perturbative) density expected from the naive Regge theory (it is reasonable 
at low with the perturbatively resummed results of a^ln^x series. 

On the other hand, the recent numerical analysis which has been done by Bliimlein 
and Vogt [|^] shows that there are no significant contributions to the evolution of gi 
from the resummation of the leading logarithmic (LL) corrections at the HERA kine- 
matical region ( x ~ 10“^). The controversial aspect between their numerical analysis 
and the assertion by Bartels et al. might be coming from the fact that the resummed 
part of the “coefficient function” is considered in Ref. [TO but not in Ref. . Bliimlein 


and Vogt did not include the resummed part of the coefficient function because this 
part turns out to fall in the next-to-leading logarithmic (NLL) corrections and de¬ 
pends on the factorization scheme adopted. It is also to be noted that the evolution, 
in general, strongly depends on the input parton densities. If one chooses a steep in¬ 
put function, the perturbative contribution will be completely washed away since the 
structure function is given as the convolution integral of the parturbative part and the 
input density. So it will be interesting to see the sensitivity of the results to the choice 
of the input densities. 

In the present paper, we numerically reanalyse the flavor non-singlet part of gi by 
taking into account the Inx resummation. We consider three different input densities: 
one is a flat density corresponding to the naive Regge prediction and others are steep 
ones in the small x region. The coefficient function can not be included consistently 
at present since the anomalous dimension has been calculated only at the LL order. 
However we consider also the effects of the coefficient function. The reason is because 
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we could firstly clarify the above controversial aspect and secondly get some idea about 
the magnitude of the NLL order corrections in the resummation approach. 

This paper is organized as follows. In section 2, we make a brief review on the 
resummation of In a; series and present an explicit expression for gi^. In section 3, 
we show our numerical results and discuss the effects of the NLL corrections. The 
interpretation of the numerical results and summary will be given in section 4. 

2 Resummation of In x terms 

The flavor non-singlet part of the polarized structure function is given by the 
formula, 

9i^{Q^,x) = ^ ( —C^^{as{Q‘^),x/y)Aq^^{Q^,y) , ( 1 ) 

z Jx y 

where is the flavor non-singlet combination of the polarized parton densities, 

2 / 2 \ 

x)=^ x) - Agi(Q^ x)) , 

and is the coefficient function, n/ is the number of active flavors with electric 
charge e^, (e^) = The perturbative evolution of the parton density is con¬ 

trolled by the DGLAP equation, 

= f —P{as{Q^),x/y)Aq{Q‘^,y) . (2) 

oQ^ Jx y 

In the above equation and in the following, we suppress the superscript NS which 
means the flavor non-singlet part. The coefficient function C{as-,y) and the splitting 
function P{as, y) are both calculable in the QCD perturbation theory. When x is hnite, 
it may be enough to compute them to the hxed-order of perturbation. In the small x 
region, however, the hxed-order calculation becomes questionable since there appear 
In” a; corrections in the higher orders of the strong coupling constant If these In” a; 
terms compensate the smallness of we must resum the perturbative series to the 
all orders to get a reliable prediction. 
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To see what terms show up at small x, it will be convenient to take the Mellin 


transform of Eq. (|^ . 


gi{Q‘^,N) = [ dxx^ ^gi{Q‘^,x) 

Jo 

= ^C(aM),N)Aq(Q\N) , 


where 


C{as{Q‘^),N) = [ x^ ^C{as{Q‘^),x) , 

Jo 

Aq{Q‘^,N) = [ x’^~^Aq{Q‘^,x) . 

Jo 

The DGLAP evolution equation Eq.(|^) becomes, 

Q^-AAq(Q\ N) = - 7 (a, W"), N)Aq{Q\ N) . (3) 

Here the anomalous dimension 7 is the moment of the splitting function, 

'y{N,as{Q‘^)) = - [ dxx^~^P{as{Q‘^),x) . 

Jo 


Eq.(|]) is easily solved to give. 


Ag(Q^ N) = Aq{Ql, N) exp 


fOs(Q^) da 


lasiQl) 


■-f{a,N) 




where (3 is the beta function. 


(3{as) = 


da.. 


(y.Q 


ainQ2 

The hrst two coefficients of the (3 function are, 

4„ „ 34 




20 


Po — “ ”3”^^ — -^CATRUf - ACpTnUf , 


with Cp = {N^ — l)/2iVc and Ca = Nc for the SU(iVc) color group and Tr = 1/2. 
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The coefficient function C{as,N) and the anomalous dimension 7 ( 0 ^ 5 , iV) may be 


expanded in the powers of 

C30 

C{a,,N) = l + ^c'^(iV)d^ 


k=l 

00 

7{a,,N) = 

k=l 


where (and in the following) we use the abbreviation, 

_ CHg 


Of o - 


47r 


The singular behaviors of the coefficient and splitting functions as x —> 0 appear as the 
pole singularities at iV = 0 in the moment space. The explicit next-to-leading order 


(NLO) calculations of the coefficient function |]T^ and the anomalous dimension |T3 
in the MS scheme show a strong singularity at iV = 0, 


7=(A') = 4(37 - + O (^) . 


( 4 ) 


whereas the leading order anomalous dimension looks like, 

^\N) = -2CFj^-CF + OiN) , 


at small N. These strong singularities (double logarithmic corrections) will persist to 
all orders of perturbative series. Indeed, at the k-th loop, the anomalous dimension 
and the coefficient function are expected to behave as, 

7(a,)^a,(^)‘ , 7a,) (5) 

Since the singularities in correspond to the ln"*“^ ^4^ singularities, our task is to 
resum these terms to all-orders in the perturbative expansion. 

Before discussing the resummed results, it may be worth mentioning the difference 
between the polarized (unpolarized flavor non-singlet) and the unpolarized flavor singlet 
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structure functions |]^. Naively one expects that the anomalous dimension behaves 
like 7 ~ at the fc-th loopQbecause there exist extra infrared and collinear 

singularities. In the case of the unpolarized flavor singlet structure functions, however, 
many of them are canceled and the true behavior at the fc-th loop is 7 ~ {ag/N)^. 
These terms can be resummed by the BFKL equation. On the other hand, above 
strong singularities survive in the polarized structure function. This fact suggests that 
the polarized structure function will receive large perturbative corrections at small x. 

The resummation of In a: singularities for the gi structure function has been done 
in Refs. 0 |^. The result for the “parton (quark) ” target with the hxed coupling 


const antp] reads. 


9i 


parton / 


.2 f-c+ioo ( Q 






N 


2 Jc-ioc 2m' J N-fo{N)/8n^' 

where /i is an arbitrary mass scale which regularizes the infrared and/or mass singu¬ 
larities. From this expression we could identify the resummed anomalous dimension 7 
and the coefficient function C to be. 


j{ag,N) = hm 7 (a„ AT) =-/o (iV)/ 87 r 2 , 

N 

C(as,N) = lim C(q!s, A^) =-- 

^ ^ ^ ^ iV-/P(iV)/87r2 


( 6 ) 

( 7 ) 


Here /q , which corresponds to the odd-signature quark-quark scattering amplitude in 
the color singlet channel, satishes the equation. 

The even-signature quark-quark scattering amplitude fg with the color octet quantum 
number in the t-channel is the solution of the equation. 


f+(N) = -871^ — — + 2Nc——f^{N) + ^^(/s+(Ar))2 , 
J8K ) ^ 87r2Ar^'^® ^ ’ 


*We follow the convention of Ref. [Q for the moment. 

Tn the genuine LL approximation, the strong coupling constant should be taken as a fixed 
parameter. 
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and given by, 


fsiN) = lQn^Ncas-^\n{e^^^^D_i/ 2 Ni{z)) with 2 ; = 


Dp{z) is the parabolic cylinder fnnction JT^. Finally we reach, 


(y.c 




1 - 


2 ir 2 ]V 


ft\N) 


Now it will be instrnctive to re-expand Eqs. (y,|3) in terms of as to see whether these 
formnlae snm np the most singnlar terms of the pertnrbative series. The expressions 
expanded np to 0{al) read, 

2 


7 = -iV 




+ 16Cf ( Cl + 2^ 


Cf 1 




+ 16Cf (^6CJ. + 12^ + 2^ - 4Cp + i + 6N, 
+ ■■■ 


cts Y 


= 

k=l 


(y.a 


ATS 

(y.g N 


( 8 ) 


C = 1 + 

+ 8 Cf (^5C^ + 8 ^ - i - 2^ 

+ 32Cf + 15^ + 2^ - 4Cp + ^^3 , ^ 

+ • ■ ■ 




as y 


k=0 


mj ■ 


(9) 


These resnlts coincide with the previons expectation of Eq.(§). Enrthermore, noting 
the relation, 

2Ca — 3Cf ~ ’ 
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which holds in SU(iVc), we can see that the resnmmed expressions Eqs.(|],^ reprodnce 
the known NLO resnlts Eqs.(^) in the MS scheme. Therefore, it is qnite plansible that 
Eqs-dlJ^) correctly snm np the “leading” singnlarities to all orders. 

Here a comment is in order concerning the scheme dependence. It is well-known 
that the anomalons dimension and the coefficient fnnction individnally depend on the 
factorization scheme and only an appropriate combination of them becomes scheme 
independent. When one considers the higher order corrections in the pertnrbation 
theory, therefore, one mnst specify the scheme adopted. Unfortnnately we do not have 
by now any appropriate factorization theorems to the problem discnssed in this paper. 
This means that we mnst be carefnl when considering the resnmmed quantities. In 
particular, the resnmmed “coefficient function” does not have any physical meaning 
until the scheme dependent part of the anomalous dimension is calculated in the same 
scheme. To clarify this issue, it is convenient to write the above results in the form 
which corresponds to the so-called DIS scheme jl^ . The DIS scheme is defined so 
that the naive parton model relation is true to all orders in perturbation theory. The 
polarized parton densities become physical observables in this scheme. The parton 
densities and anomalous dimension in the DIS scheme are obtained by making the 
transformations. 


Aq Aq^^^ = CAq , 

C/O^s 

Using the resnmmed 7 and C Eqs.(H,H), we get the resnmmed part of the anomalous 
dimension in the DIS scheme. 
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DIS 


= A'f 7 -= f V + AA'Wf 4) ‘ + o f ivhll) 


k=l 


k=2 


Ny 


NV 


( 10 ) 


where the second terms come from the resnmmed coefficient function and are nu¬ 
merical numbers independent of N. The above equation tells us that the resnmmed 





coefficient function belongs to the NLL order^ corrections in the context of the resum¬ 
mation approach. Then, one must include the NLL order anomalous dimension which 
has not yet been available to see the effects of the coefficient function. This is the 


reason why the authors in Ref. |11| throw away the coefficient function. 


3 Numerical Analysis 

Numerical analysis of the spin structure function in the small x region was done in 
the context of the small x resummation approach in Ref. m- They obtained the result 
that the small x resummation effect is not signihcant despite of a naive expectation 
discussed in Ref. |T^. In this section, we numerically reanalyze the behavior of 
structure function to show how the hnal results are sensitive to the choice of the input 
parton densities. In conjunction with the claim in Ref. |^, we also consider the effects 
from the resummed coefficient function. As already discussed in section 2 , we can not 
include the coefficient function in a theoretically consistent way. However we believe 
that the inclusion of the coefficient function could shed some light on the size of the 
NLL corrections in the resummation approach. 

At hrst, we explain our method to estimate the g^^ structure function numerically. 
Our starting point is the expression, 

NS(i^2 ^ rc+ico dN j^ / /-"hQ ) das njs\ nS(i^2 at', 

9i {Q,x)= —X exp - / —7 g^ [Qo^^) ■ 

J C—iOQ Y j Ois(Q^) P J 

The anomalous dimension which includes the resummation of In”a; terms is orga¬ 
nized as follows, 

y^^^{N)=as'y\N) + aly^{N) + K{N,as)-/3-^\n{l + asC^ + H{N,as)) , ( 12 ) 

where 7 ^’^ and are respectively the usual anomalous dimension and coefficient func¬ 
tion at the one and two-loop hxed order perturbation theory. K{N,as) {H{N,as)) is 
^This fact implies the LL resummed anomalous dimension 7 being scheme independent. 
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the resummed anomalous dimension Eq.(|^) (Eq.(^) with fc = 1, 2 (/c = 0,1) terms be¬ 
ing subtracted because those terms have already been included in the usual anomalous 
dimension and coefficient function. 

K{N,as) ^ 

It should be noted here that the anomalous dimension at = 1 plays a special role 
for the non-singlet gi structure function. In a language of the operator product expan¬ 
sion, 'j{N = 1) is the anomalous dimension of the (non-singlet) axial vector current. 
Since the (non-singlet) axial vector current is conserved, the corresponding anomalous 
dimension should vanish. The perturbation theory guarantees this symmetry order by 
order in the expansion. However, the resummation of the leading singularities in 
N does not respect this symmetry. Therefore, we need to restore this symmetry “by 
hand”. In this paper, we multiply K{N, a*) by (1 — N) Q. 

K(N,a,) ^ K(N,a.)(l - N) . 


Of course, this is not a unique prescription and one can choose other procedure^ which 
satisfies the condition of limAr^i K{N, ag) = 0. 

Now let us explain how to perform the Mellin inversion Eq. ([Tl|) which is the integral 

in the complex iV-plane. At first, we must know the Mellin transform of the input 

function gi{QQ,N). It is easy to obtain an analytical form for it in the complex N- 

plane since we assume a simple function (see below) for the input density. Next we 

need an analytically continued expression of the anomalous dimension in the 

complex A^-plane. For the gi structure function, only odd moments are defined. So 

we replace (—1)^ by (—1) in the expression of the anomalous dimension obtained in 

Ref. 0. The integration contour in the Mellin inversion should be on the right of the 
^Our final conclusion remains the same qualitatively if we choose other prescription. 
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rightmost singularity of the integrand. The contour integration along the imaginary 
axis from c — ioo to c + too is numerically inconvenient due to the slow convergence 
of the integral in the large |A^| region. To get rid of this problem, we deformed the 
contour to the line which have an angle 0 (0 > 7r/2) from the real N axis. By this 
change of the contour, we have a damping factor exp(|iV| ln(l/x) cos^) which strongly 
suppresses the contribution from the large \N\ region. In the integration along this 
new contour, we will be able to cut the large |A^| region. Finally we have checked the 
stability of results by changing the contour parameter. One can find the details of this 


technique in Ref. [18 


We choose the starting value of the evolution to be Ql = 4GeR^. We calculate the 
evolution for three types of the input densities A, B and C: A is a function which 
is flat at small x {x°‘, a ~ 0), B is slightly steepp| (a ~ —0.2) and C rises more steeply 
(a ~ —0.7). The explicit parametrization used in this paper is l]^, 

Aq(Ql, x) = 77(a, /3, a)rix°‘{l — xY{l + ax) , 

where A^ is a normalization factor such that / dxNx°'{l — x)^{l + ax) = 1 and rj = 

{qa = 1.26) in accordance with the Bjorken sum rule. A, B and C correspond to the 
following values of parameters. 


A : a = +0.0 , /3 = 3.09 , a = 2.23 , 

B : a = -0.2 , f3 = 3.15 , a = 2.72 , 

C : a = -0.5 , f3 = 2.41 , a = 0.02 . 


In our analysis we put the flavor number rij = 4 and Aqcd = 0.23GeI7. 

First we estimate the case which includes only the LL correction 7 . The evolution 

kernel in this case is obtained by dropping H{N,as) in Eq.(^)|[J]. This is a consistent 

^This choice is essentially the same as one in Ref. []lH . 

II We take into account the fixed-order NLO corrections exactly. 
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approximation in the resummation approach. Fig.la (lb, Ic) shows the results (dashed 
curves) after evolving to = 10,10^, lO'^GeV^ from the A (B, C) input density (dot- 
dashed line). The solid curves are the predictions of the NLO-DGLAP evolution. 
These results show a tiny enhancement compared with the NLO-DGLAP analysis and 
are consistent with those in Ref. HQ In the case of G, we can not discriminate 
a difference between the LL and DGLAP results. The enhancement is, as expected, 
bigger when the input density is flatter. However any signihcant differences are not 
seen between the results from different input densities. 

Next, we include the NLL corrections coming from the resummed “coefficient func¬ 
tion”. We show the results in Fig.2 by the dashed curves. (Other curves are the same 
as in Fig.l.) The results are rather surprising. The inclusion of the coefficient function 
leads to a strong suppression on the evolution of the structure function at small x. 
Since the effects from the coefficient function fall in the NLL level, the LL terms are 
expected to (should) dominate in the small x. However our results imply that the 
LL approximation is not sensible in the small x region we are interested in. As the 
resummed coefficient function is only a part of the NLL correction, we can not present 
a dehnite conclusion on the (full) NLL correction. But it is obvious that the NLL cor¬ 
rection is very important at the experimentally accessible region of x. In the following 
section, we explain why the coefficient function leads to such suppression. 

4 Discussion and Summary 

In the previous section, we have shown that although the LL resummed effect is very 
small at the experimentally accessible region of x, a part of the NLL resummed con¬ 
tribution from the coefficient function drastically changes the predictions. 

To understand these numerical results, it will be helpful to remember the per- 

**We have also calculated gi with the input function used by Bliimlein and Vogt and could reproduce 
their results. 


12 



turbative expansion of the resummed anomalons dimension and coefficient fnnction 
Eqs-dlJ^). By using the explicit values Nc = 3, Cp = 4/3, we obtain for the anomalous 
dimension in the DIS scheme Eq.(^), 


7 


DIS 


= N 


- 0.212 


as 


0.068 


as 

Ip 


0.017 I Jl) - 0.029 (||) + 


+ 

+ ■■■ 


0,141 I +0,119 (^)’ +0,069 + 


(13) 


Here note that: (1) the perturbative coefficients of the LL terms (the hrst part of 
Eq.(^)) are negative and those of the higher orders are rather small number. This im¬ 
plies that the LL corrections push up the structure function compared to the hxed-order 
DGLAP evolution, but the deviations are expected to be small. (2) the perturbative 
ones from the NLL terms (the second part of Eq. (P^ ), however, are positive and some¬ 
how large compared with those of the LL terms. This positivity of the NLL terms 
has the effect of decreasing the structure function. This fact that the coefficients with 
both sign appear in the anomalous dimension should be contrasted with the case of 


the unpolarized structure function |1T9 |. 

Now it might be also helpful to assume that the saddle-point dominates the Mellin 
inversion Eq.([TI|). We have numerically estimated the approximate position of the 
saddle-point and found that the saddle-point stays around iVsp ~ 0.31 in the region of 
X ~ 10“® to 10“^. (Of course the precise value of the saddle-point depends on x, Ql 
and Q^.) By looking at the explicit values of the coefficients in Eq. (p!3D , the position 
of the saddle-point seems to suggest that the NLL terms can not be neglected. Since 
the coefficients from the higher order terms are not so large numerically, it is also 
expected that the terms which lead to sizable effects on the evolution may be only 
hrst few terms in the perturbative series in the region of x we are interested in. We 
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have checked that the inclusion of the hrst few terms in Eq. ([T3| ) already reproduces the 
results of section 3. Fig.Sa (3b) shows the numerical results of the contribution from 
each terms of the NLL corrections in Eq.(|T^) at with the A (B) type 

input density. The solid (dot-dashed) line corresponds to the NLL (LL) result. The 
long-dashed, dashed and dotted lines correspond respectively to the case in which the 
terms up to the order a^, are kept in the NLL contributions. One can see that 
the dotted line already coincides with the full NLL (solid) line. These considerations 
could help us to understand why the NLL corrections turns out to give large effects on 
the evolution of the gi structure function. 

The hnal discussion concerns the convergence issue of the perturbative series. As 
discussed in Refs. [^, one must be careful when applying the perturbative ap¬ 
proach to small X evolution. The integrand in Eq. (pT]) has a singularity in the moment 
space. This (rightmost) singularity is equal to that of /o~(A^). The numerical value A^^o 
of the singularity position is A^^o ~ 0.304. This means that the N can not become so 
small. On the other hand, the approximation scheme in the resummation approach is 
sensible only for small N. This apparent contradiction will be solved by analyzing the 


evolution in x space By explicitly solving the evolution in x space, it has been 
pointed out that the saddle-point method is not a good approximation in the case 
of the unpolarized structure function. Although we have not used the saddle-point 
approximation to solve the evolution, the previous explanation relying on this method 
can be misleading. So according to Refs. ||^ |^, we have tried to solve the evolution 
in X space with hrst several terms of the perturbative expansion being kept and what 
we found is that the conclusion does not change. The numerical results are essentially 
the same as Fig.3. 

In summary, we have performed numerical studies for the havor non-singlet gi 
structure function at small x by incorporating the all-order resummed anomalous di- 
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mension and a part of the NLL corrections from the resnmmed coefficient fnnction. 
Onr results show that the resummed coefficient function has an effect which suppresses 
the structure function at small x. Including only the resummed coefficient part is not 
theoretically consistent, and so one should take into account also the anomalous dimen¬ 
sion at the NLL level. However, our results suggest that the LL analysis is unstable, in 
the sense that a large suppression effect comes from the resummed coefficient function 
which should be NLL correction. We have explained why the inclusion of a part of 
the NLL corrections leads to such unexpected results. We need a full NLL analysis 
to make a definite conclusion on whether the all-order resummation approach predicts 
a rise of the flavor non-singlet gi structure function in the experimentally accessible 
small X region. 
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Figure Captions 


Fig. 1 

The LL evolution as compared to the DGLAP results with the flat input A (Fig. la) 
and steep ones B (Fig. lb) and C (Fig. Ic). 

Fig. 2 

The NLL evolution as compared to the DGLAP results with the flat input A (Fig. 2a) 
and steep ones B (Fig. 2b) and G (Fig. 2c). 

Fig. 3 

Gontributions from the fixed order terms in the NLL resummation with the flat input 
A (Fig. 3a) and steep one B (Fig. 3b). 
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